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^invention concerne un procede de localisation d'un ou de 
plusieurs emetteurs. 

II s'applique notamment a la localisation d'emetteurs au sol a partir 
d'un engin mobile sans avoir de connaissance a priori sur les signaux emis, 

5 Le domaine technique est notamment celui de la localisation 

passive d'emetteurs. 

La figure 1 illustre une localisation aeroportee. L'emetteur est a la 
position (x 0 ,y 0 ,zo). Le porteur a i'instant t k est a la position (x k ,y k ,z k ) et voit 
l'emetteur sous incidence (6(t k ,Xo t yo,zo), A{t k ,x 0 ,yo,z 0 )). Les angles 
10 6(t,xo,yo>zo) et A(t,x 0 ,yo,z6) evoluent au cours du temps et dependent de la 
position de l'emetteur ainsi que de la trajectoire du porteur. 

Les angles tffrxa,yo,Zo) et A(t,x 0 ,y 0 ,Zo) sont reperes par rapport a 
un reseau de N antennes pouvant etre fixe sous le porteur comme le montre 
la figure 2. 

15 II existe actuellement plusieurs techniques permettant |de 

determiner la position (x^, Ym. z m) d ' un emetteur. Ces techniques. de 
localisation different notamment par les parametres qui sont estimes en 
instantanee au niveau du reseau des capteurs. Elles peuvent etre classees 
de la maniere suivante: 

20 Utilisation de la goniometrie, 

Ces techniques sont connues et utilisees dans Tart anterieur;- 
Dans la plupart des cas, elles sont basees sur une goniometrie ID en 
azimut Les azimuts e^^x^y^Zm) associes au ni*™ emetteur sont 
mesures pour differents instants t k . En utilisant la position (x kj yk,z k ) du porteur 
25 a Tinstant correspondant k, une position {x mk ,y mk ,z mk ) de I'emetteur m est 
estimee par une intersection sol. La position (Xk,y k ,z k ) du porteur est donnee 
par un GPS, son orientation est obtenue par un compas dans le cas d'un 
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porteur terrestre et par une centrale de navigation dans le cas d'un aeronef. 
A partir de toutes les positions (x mk ,y mk ,z mk ), la methode effectue une 
extraction de donnees permettant de determiner les M positions dominantes 
(x m ,y m ,z m ) des emetteurs incidents. La localisation est obtenue par 
triangulation ou par intersection sol (goniometrie 2D). L'inconvenient des 
techniques de triangulation est qu'elles necessitent un defilement important. 
D'autre part, les techniques de goniometrie doivent utiliser un reseau de 
capteurs non ambigue pour fournir les incidences. Ceci a pour inconvenient 
de necessiter une table de calibration et de limiter la taille du reseau de 
capteurs et par consequent de fournir des incidences limitees en precision. 

Utilisation de la difference de phase entre 2 capteurs eloignes, 

La difference de phase ^f fc x 0) y 0 ,z 0 ) entre capteurs depend de la 
positions des 2 capteurs ainsi que de I'incidence (^ fe x 0> y 0|2o ), J(f fo x 0 ,y 0> z 0 )) 
de I'emetteur. Cette phase qui depend du temps est directement liee a la 
position (xo.yo.zo) de I'emetteur. En consequence, en etudiant la fonction du 
temps A<p(t,Xo,yo,zo) il est possible d'en d<§duire la position (x 0l y 0 ,z 0 ) de 
I'emetteur. Dans cette famille d'application les 2 capteurs sont eloignes pour 
augmenter la precision de la mesure de la phase. Ceci a pour inconvenient 
de faire varier la difference de phase A<p(t,x 0 ,y 0 ,z 0 ) en fonction du temps sur 
plus de 2n et la technique necessite alors une etape permettant de derouler 
la phase sur plus de 2tc D'autre part dans cette technique la phase est 
mesuree en effectuant directement une intercorrelation entre 2 capteurs ce 
qui ne permet pas de traiter le cas multi-emetteurs. 

Utilisation de la mesure de la frequence porteuse de I'emetteur, 

Ces techniques exploitent le fait que la frequence porteuse 
estimee est la somme de la frequence porteuse de I'emetteur et du decalage 
doppler du a la Vitesse de deplacement du porteur. Le decalage doppler a 
I'avantage de dependre de la position (x 0 ,y 0 ,z 0 ) de I'emetteur et d'etre aussi 
une fonction du temps Af(/,x 0 ,y 0 ,z 0 ). En consequence en etudiant la fonction 
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du temps Af(f,X0,y o ,z o ) il est possible d'en deduire la position (x 0 ,yo,za) de 
I'emetteur. La mesure de ce decalage doppler presente toutefois comme 
inconvenient de necessiter des emetteurs ayant des formes d'ondes 
particulieres. Cette mesure de frequence peut se faire par des techniques 
5 cycliques supposant que le signal emis est non circulaire. 

Utilisation des temps de propagation, 

Ces techniques exploitent les differences de temps de propagation 
entre aeriens (TDOA ou Time difference of arrival) qui sont directement liees 
aux distances respectives de I'emetteur aux differents aeriens et done a la 

10 position (x 0 ,yo,zo) de I'emetteur. En utilisant au moins trois aeriens 
suffisamment espaces, il est possible de deduire la position (x 0t yo,Zo) de 
I'emetteur par localisation hyperbolique. ^inconvenient de ces techniques est 
qu'elles ne peuvent etre mises en ceuvre en contexte mono porteur en raison 
des espacements considerables requis entre aeriens. D'autre part dans ces 

15 techniques, la difference de temps est mesuree en effectuant directement 
une intercorr^lation entre 2 capteurs, ce qui ne permet pas de traiter le cas 
multi-6metteurs. 

Le procede selon I'invention repose notamment sur une 
nouvelle approche d'estimation directe des positions (x m , y m , z m ) de chacun 
20 des emetteurs a partir d'une analyse parametrique du signal multi-voies a 
divers instants t k sur une duree Dt. Uanalyse parametrique a notamment 
pour fonction supplementaire de separer les differents emetteurs a chaque 
instant t k . On associe ensuite les parametres d'un meme emetteur issus des 
differents instants t k pour finalement localiser chacun des emetteurs. 

25 LMnvention concerne un procede de localisation d'une ou de 

plusieurs sources, la ou lesdites sources etant en mouvement par rapport 
a un reseau de capteurs, le procede comportant une etape de separation 
des sources afin d'identifier les vecteurs directeurs associes a la reponse 
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des capteurs a une source d'incidence donnee. II est caracterise en ce 
qu'il comporte au moins les etapes suivantes : 

o associer les vecteurs directeurs a lm ...a Km du m'eme emetteur obtenus 
respectivement aux instants t,...^, 

• localiserle m'eme emetteur a partir des vecteurs a 1m ...a,, associes. 

Le precede selon I'invention presente notamment les avantages 

suivants : 

• il permet de localiser en plus de la position en (x, y, z) d'un emetteur son 
vecteur vitesse, 

o il s'applique lorsque Ton est en presence de un ou plusieurs emetteurs 
incidents, 

o sa mise en oeuvre ne necessite pas de connaissances particulieres sur le 
signal emis, 

» il permet d'utiliser un reseau de capteurs ambigue (c'est-a-dire tel que 
plusieurs incidences peuvent etre associees a la meme reponse du 
reseau) qui ont I'avantage d'etre grand et ainsi d'etre plus precis et plus 
robuste aux phenomenes de couplage entre aeriens ou plus 
g§neralement aux erreurs de mod^lisation du reseau d'aeriens, 

• il peut dtre mis en oeuvre sur des reseaux calibres en (9,A). 

• il peut etre mis en oeuvre sur des reseaux a antennes a diversites 
d'amplitude comme les antennes colocalisees : reseau avec des dipoles 
de meme centre de phase et ayant des orientations differentes. 

D'autres caracteristiques et avantages de I'objet de la presente 
invention apparartront mieux a la lecture de la description qui suit donnee a 
titre illustratif et nullement limitatif a la lecture des figures annexees qui 
repr^sentent : 
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• la figure 1 le schema de principe de la localisation d'un emetteur de 
position au sol au moyen d'un aeronef, 

» la figure 2 la relation entre un reseau d'antennes et I'incidence d'un 
§metteur, 

5 © la figure 3 un schema general expliquant le fonctionnement du procede 
selon I'invention, 

o les figures 4, 5 et 6 des exemples de mise en oeuvre du procede selon 
(invention. 

Afin de mieux faire comprendre I'objet de la presente invention, la 
10 description qui suit est donnee a titre illustratif et nullement limitatif pour 
localiser plusieurs emetteurs disposes au sol au moyen d'un reseau de 
capteurs equipant un aeronef en mouvement. Un tel systeme est par 
exemple decrit a la figure 1 . L'aeronef est equipe d'un processeur adapte a 
mettre en ceuvre les etapes du procede selon I'invention. 

15 Le procede peut aussi etre mis en oeuvre dans le cadre de 

vehicule en mouvement au sol. 

La figure 3 represente, dans un diagramme temps-amplitude du 
signal, le signal x(t) compose d'une combinaison des signaux des emetteurs 
a differents instants t^ t 2 , ...tK- Sur cette figure les differentes etapes mises 
20 en ceuvre a savoir la separation des emetteurs SE et Testimation 
parametrique EP, Tassociation des parametres de chaque emetteur, la 
localisation d'un Emetteur sont resumees. 

En presence de M emetteurs, le procede dispose, a Tinstant t en 
sortie des N capteurs du reseau, du vecteur x(f) representatif du melange 
25 des signaux des M emetteurs. Autour de I'instant k le vecteur x(t+t k ) de 
dimension A6c1, representant le melange des signaux des M emetteurs, 
s'exprime de la maniere suivante : 
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*(t+tk) =g a(Q km , A km ) s m (t+t k )+ b(M k )=A k s{t+t k ) + b{t+t k ) pour \t\<At/2 (1) 



ou b(0 est le vecteur bruit suppose gaussien, a(6,A) est la reponse du reseau 
de capteurs a une source d'incidence (e,A), A^ap*. A„)... a^y, 
s(f)=[si(r)-.s iW (0] T , Qim=6{t kt x m ,y ms z m ) et ^fcm=A{t k ,x m ,y m ,z m ). Dans ce modele, 
la matrlce A k de melange depend de I'instant t k d'observation. 

5 Le vecteur directeur de I'incidence correspondant au rri° me emetteur k 
I'instant t k 

akm= aiQkm, A km )=a(t k ,x m ,y m ,z m ) du ni hme Emetteur (2) 

est une fonction connue de t k et de la position de I'emetteur (x m ,y m ,z m ). 

Le procede selon I'invention comporte au moins les etapes 

suivantes : 

10 1. estimer un ou plusieurs parametres associes a la position de la source, 
par exemple les vecteurs directeurs, les incidences, la position, etc. et 
separer les M emetteurs pour les differents instants t k , ce qui consiste a 
identifier les vecteurs directeurs d'incidence a*™ pour (1<m</W). Cette 
premiere etape est par exemple effectuee par des techniques de 

15 separation de sources connues de I'Homme du metier, 

2. associer les parametres estimes pour le m'eme emetteur, par exemple en 

associant les differents vecteurs directeurs d'incidences, a 1m a km 

obtenus respectivement aux instants t-\ t ..\, 

3. localiser le rri*™ emetteur a partir des vecteurs associes. 
20 Etape dissociation 

En presence de M emetteurs et apres separation de sources, le 
procede possede a I'instant t k les M signatures a^ pour (1<m<M). A I'instant 
t M la separation de source donne les M vecteurs b/ pour (1</</W). L'objectif 
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de ce suivi est de determiner pour le rri* me emetteur, I'indice i(m) qui 
minimise I'ecart entre a/on et h,-^;- Dans ce cas on en deduira que a k + 1r ^= 
b i(m) . Pour effectuer cette association on definit par exemple ia distance entre 
deux vecteurs u et v par : 



u H v\ 2 

d(u. v) =1- 



(3) 



(u H u\v H v) 

5 Ou u H est le transpose conjugue du vecteur u. 
Dans ces conditions I'indice i(m) verifie : 
d(a fcm , h, (m} ) = min [d(a km , &,)] (4) 

done 



H I 2 
\ a hn b i(m)\ 



d(a to , b l(m) ) = mm[1- r H V u b y 



Dans cette association on considere une fonction fi m associee au ni eme 
10 emetteur: 

Mfc) = dCaton, a 0m ) (5) 

Au fit de I'association on obtient pour chaque emetteur m et pour 1<m<M, la 
fonction p m (f). Cette fonction a notamment pour objectif d'eliminer les instants 
t k dont la valeur 0 m (f*) paraTt trap eloignee d'une interpolation de la fonction 
MO. e'est-a-dire que I'on elimine les instants aberrants qui peuvent etre 

15 associes a d'autres emetteurs. On definit une zone de tolerance +/-A autour 
de la courbe definie par la fonction Mfc). Cette zone de tolerance dependra 
de la precision d'estimation des vecteurs directeurs a*™. En particulier en 
presence de M=1 source la zone sera de I'ordre de A=3/V5At (oil At est le 
temps elementaire d'estimation parametrique illustre fig.3 et B est la bande 

20 instantanee du signal x(t) ). 
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Les etapes de cette association pour K instants t k sont par 
exemple les suivantes : 

Etape ASE - 1 : Initialisation du processus a k=2. Le nombre M d'emetteurs 
initial est par exemple determine par un test de detection du nombre de 
5 sources a I'instant t 0 connu de I'Homme du metier, 

Etape ASE - 2 : Pour 1<m<M determination des indices i(m) en appliquant 
I'equation (4) et en utilisant le vecteur a Km avec 1 <m<M et les vecteurs b, 
identifies a I'instant pour (1</</W), 

Etape ASE - 3 : Pour 1 <m<M effectuer ('operation a k+1 ^ b l(m)i 

10 EtapeASE-4: Incrementation /«-/c+1 et si /c<Kretour a I'etape ASE-1, 

Etape ASE - 5 : A partir de la famille d'instants <£={ u<...< t K ), eliminer les / 
instants fc® tel que les coefficients 0 m (Q n'appartiennent pas a la zone 
delimitee par la courbe ^interpolation des j3 m (t k ) et la zone de tolerance A. 
On eliminera aussi les instants t* ou \U%)- Mt k -,)|<A. Apres ce tri la 
15 nouvelle famille d'instants est 0={ *<...< t, } et on pose /C=/. 

A la fin de ces etapes, le precede a determine les vecteurs a 1m 
a/c^associes au ni bme emetteur. 

Localisation d'un emetteur 

Le precede determine la position du rri°™ emetteur a partir des 
20 composantes des vecteurs a 1m jusqu'a Ces vecteurs a^ ont la 
particular^ de dependre de I'instant t k et surtout de la position (x^z,,,) de 
I'emetteur. En particulier pour un reseau compose de N=2 capteurs espaces 
d'une distance de dans I'axe du porteur le vecteur verifie a km : 



3km = 



, 1 



(6) 

= a(^x m ,y m ,z m ) 
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La valeur 1 de la premiere composante correspond au capteur de 
reference. D'apres la figure 1, I'incidence (0(t k ,x m ,y m ,z m ), A(t k ,x m ,y m ,z m ) ) peut 
etre directement calculee a partir de la position (x k ,y k ,z k ) du porteur a I'instant 
t k et la position (x m ,y m ,z m ) de I'emetteur. 

5 Etape de transformation du vecteur 

Selon un premiere variante de realisation, le procede comporte 
une etape de correction des aim, La mesure des vecteurs directeurs a km est 
genSralement obtenue a un facteur complexe pres indetermine. Selon cette 
premiere variante, le procede comporte une £tape qui consiste a changer la 
10 reference de phase du vecteur directeur mesure en se ramenant au 
ban/centre de phase (defini a un coefficient scalaire constant pres que Ton 
peut fixer a 1). Cette operation, est realisee, par exemple, en estimant le 
coefficient de correction determine par la transformation suivante des a fcm 
en a' Km : 



i 




Le coefficient de correction n'est pas totalement determine par 
cette expression compte tenu de I'indStermination d'ordre N de la racine 
complexe. Un suivi de Involution de phase pendant la p<§riode d'observation 
20 est done effectue. 

Le coefficient complexe etant d<§fini a un facteur pres parmi les N 
racines N l6me de I'unite, le suivi de phase consiste a fixer arbitrairement le 
premier coefficient de correction (en prenant la racine 1 par exemple), puis a 
determiner a chaque nouvelle iteration k+1, le coefficient qui minimise les 
25 ecarts de phases moyens entre le vecteur directeur recentre a k+1 et le 
vecteur recentre a I'instant k. 
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Le critere de minimisation, pour des mesures a la meme 
frequence, peut §tre egal a : 



™B Zmin(mod( 



arg( 



«. (0 



ZnXIn - mod( 



arg( 



«.(0 



(8) 



ou les a k+1 sont les vecteurs directeur recentres avec le coefficient de 
correction determine arbitrament par I'une quelconque des racines Nemes 
de Pexpression. Pour des mesures a des frequences differences, il est 
possible de comparer les phases des composantes des deux vecteurs 
directeurs en les corrigeant d'une puissance donnee par le rapport de ces 
deux frequences. 



Si I'on considere les vecteurs b km = a' km , il est alors possible de 
comparer cette mesure a la valeur theorique b(t k ,x m ,ym,z m ) pour laquelle le 
vecteur directeur theorique a(t kl x m ,y m , 2m ) est calcule pour une origine 
consideree au barycentre (geometrique) theorique de phase (lieu 
geometrique pour lequel la somme theorique des differences de phase 
15 s'annule). Ce lieu ne coincide pas (en general) avec le centre de phase du 
reseau. 



Selon une autre variante de realisation, le procede comporte une 
etape de transformation du vecteur a to en un vecteur dont les 
composantes sont formees a partir des composantes du vecteur a**,. En 
particulier, le procede construit par exemple le vecteur de dimension (N- 
1)x1 en choisissant un capteur de reference en n=i : 
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»ta(i)/«*.(0 

»i»('-l)/a*.(0 



= b(f to x m ,y m ,z m ) 



(9) 



ou aton(f) est la / 6me composante de 

Les composantes de b* m correspondent dans ce cas aux rapports des 
composantes du vecteur a^n et du vecteur a^i). 

Ainsi dans I'exemple de I'equation (6) en fixant *=1 on obtient : 



bfcm = 



exp^72^^cos(^ , x m ,y m ,z m ))cos(A{t k ,x m , y„„z m ))j (10) 



= a km (2)/ a fcm (1) 

Sachant que les vecteurs directeurs a*m sont estimes avec une 
5 certaine erreur tel que a*^ a(k,x mi y m ,z m )+ e^, on peut en deduire qu'il 
en est de meme pour le vecteur transforme h km de (9) . 

Etaoe de maximisation d'un critere de correlation 

Sachant que le vecteur a*m est une fonction de la position 
(Xm.ym.Zm) de I'emetteur il en est de meme pour le vecteur b^- Le procede 
10 comporte une 6tape de maximisation d'un critere de correlation vectorielle 
normalise U(x,y,z) dans I'espace (x,y,z) de position d'un emetteur ou 

\*>K Hv K(. x >y> z )\ 2 



Avec 
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(11) 



et W/c = 



w 



w 



Km 



Le vecteur bruit w* a pour matrice de covariance R=E[w K w K H ] En 
faisant I'hypothese que la matrice R est connue, ie critere peut etre envisage 
avec une technique de blanchiment. Dans ces conditions on obtient le critere 
U'(x,y,z) suivant : 



L/f'(x,y,z)= r — - — -u^ ~ - -i 

Avec R=E[w K w / c H ] 



(12) 



5 . II faut remarquer que les criteres des equations (11) et (12) sont egaux 
lorsque R=o 2 l, c'est a dire lorsque les erreurs sont considers de niveau 
egal sur tous les capteurs et independantes entre capteurs. Le critere 
Ux,y,z) de ('equation (12) est done valable pour un bruit de modele w* de 
statistiques blanches. 

10 o Les criteres L^x.y.z) et L K >(x,y,z) sont compris entre 0 et 1 et verifient 
Mx,y,z)= Uc'(x,y,z)=1 pour la position (x^aO du nF» emetteur. Cette 
normalisation permet de fixer un seuil de bonne localisation n . Ainsi tous les 
maximums (x m ,y mj z m ) de L^y.z) qui verifient L K (x mi y fni z m )> „ sont 
cons.deres comme des bonnes localisations. Le seuil peut etre fixe en 

1 5 fonction d'une connaissance approchee des statistiques de w* 

• Les criteres L^y.z) et L^(x,y,z) ont I'avantage de pouvoir mettre en 
ceuvre une technique de localisation en presence d'un reseau de capteurs 
calibres dans I'espace (e,A). Sachant qu'a I'instant t k on connaTt la relation 
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analytique liant I'incidence (e(f*,x,y,z), A(fo,x,y,z)) de I'emetteur a sa position 
(x.y.z), on peut alors deduire a partir de i'incidence (6(fetX,y,z), A^x.y.z)) le 
vecteur a(f to x m ,y m ,z m )= a(e(k,x,y,z), A^x.y.z)) en realisant une interpolation 
de la table de calibration (relative aux antennes calibrees). Remarquons 
5 cependant que ce procede est insensible a un biais en phase (en raison du 
critere de correlation vectorielle). 

• Ces criteres permettent aussi de tenir compte de la phase et de I'amplitude 
des composantes de a(9,A). La methode peut done etre envisagee avec des 
reseaux a antennes colocalisees a diversity de diagramme. 

10 II faut remarquer que dans un contexte aeroporte la connaissance de 
I'altitude h de I'avion permet de reduire le calcul du critere dans I'espace de 
recherche (x,y) en posant z=h. Dans I'exemple des equations (6) et (10) le 
vecteur v K {x,y,z) s'ecrit de la maniere suivante : 



VK(x,y,z) = 



exJ^J2xj cos(0fo , x, y, z)) cos(A(;, , x, y, z))j 
exp^y 2w | cos(0{t K , x, y, z)) cos(A(^ , x, y, z)) J 



(13) 



Dans ce procede il est possible d'envisager une etape 
15 d'initialisation de I'algorithme a K=Kq, puis ensuite de calculer de fagon 
recursive le critere U(x,y,z). Dans ces conditions LK(x,y,z) se calcule 
recursivement de la facon suivante : 



U + i(x,y,z)= 



\a K »{x,y,zf 



ou 



aK + i(x,y,z)= a^x.y.z) + b K+ i m H b(f K+ i,x,y,z) (14) 
YK+i(x,y,z)= YK(x,y,z) + b(fK + i,x,y,z) H b(fo + i,x,y,z) 
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ftai=pV + b K+ i m H b K+1 



m 



Les coefficients a K ^x,y,z)= c^x.y.z) Y K + i(x,y,z)= Y *(x,y,z), p> +1=p> 
sont des spectres intermediates permettant de calculer ^(x.y.z). 

Lorsque les vecteurs b(^i,x,y |2 ) et bton sont de normes 
constantes egale ap\a relation de recurrence de I'equation (14) devient : 

u +l( x >y ,z)= fc ? (*»< 

°* (15) 
ot/c + i(x,y,z)= a^x.y.z) + b K+1m H b(k +1 ,x,y,z) 

5 Le precede est decrit jusqu'ici en supposant que les emetteurs ont 

des positions fixes. II peut facilement s'etendre au cas de cibles mobiles de 
vecteur vitesse (v xrtl) v ym ,v zm ) pour lesquelles on dispose d'un modele 
devolution. Dans ces conditions, I'incidence du rri km6 emetteur se 
parametrise de la maniere suivante (1 6): 

Qknrm, x m - Vwn t k ,y m -v ym t k , z m - Vzm t k ) et A^Ai f to x m - Vxm t k ,y m -v ym t k , Zn7 -v zm t k ) 

10 ou (x m ,y m ,z m ) est la position de I'emetteur a I'instant t 0 et ( Vxm , Wzm ) les 
composantes de la vitesse de I'emetteur a I'instant t 0 . Dans ces conditions le 
vecteur de I'equation (9) est parametre par (x m ,y m ,z m ) et (v^.vg de 
la manidre suivante : 

bfen = b(ffoX /T , f y m) z /n( V xm ,V ym ,V 2n „)+ Wfen (17) 

De maniere naturelle les criteres de localisation L K et U' des 
15 equations (11) et (12) ne sont plus parametres seulement par (x.y.z) mais 
aussi par (v x ,v y ,v 2 ). Le procede consiste done a maximiser le critere L^x.y.z, 
v x ,v y ,v 2 ) en fonctions des 6 parametres (x,y,z, v x ,v y> v z ). 

Le procede peut s'appliquer a un tres grand nombre de mesures. 
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Dans ce cas, le procede comports une etape de reduction de la complexity 
numerique de calcul (qui est fonction du nombre de mesures) en diminuant 
K. Le procede prevoit d'effectuer sur les mesures elementaires les 
traitements suivants : 

5 - decimation des instants t k , en eliminant les instants voisins pour lesquels 
revolution de la courbe Mfe) n ' est P as significative, 

- filtrage (lissage des mesures qui sont les vecteurs directeurs) et sous 
echantillonnage, 

- les mesures sont ensuite fusionnees sur une duree definie (extraction par 
10 association de vecteur directeur pour produire une mesure de synthese). 

Recapitulation des etapes du procede 

Le procede de localisation de plusieurs emetteurs utilisant /("instants t k 
peut se resumer par les etapes suivantes : 

Etape n°1 : Identification des vecteurs pour (1<m<M) aux K instants f k en 
15 appliquant par exemple une technique de separation de sources et 
d'identification de sources comme decrit dans les references [2] [3]. 

Etape n°2 : Association des vecteurs a 1/n jusqu'a a^ obtenus aux instants 
respectifs ft... t K associes au rri*"™* emetteur pour 1</n<Men appliquant les 
etapes ASE-1 jusqu'a ASE-5 decrit ci-dessus. 

20 Etape n°3 : Initialisation du processus a m='\ 

Etape n°4 : Transformation des K vecteurs a*™ en des vecteurs b Km comme 
le suggere I'equation (9). 

Etape n°5 : Calcul et maximisation du critere L^x.y.z) de I'equation (1 1) pour 
obtenir la position (x m ,Vm,Zm) du ni bme emetteur. 

25 Etape n°6 : Incrementation m<-m+1 et si m<M retour a I'etape n°3 
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Afin d'affiner I'estimation de la position ( Wm , ym ) des emetteurs les etapes 
du procede peuvent etre realisees de facon iterative de la maniere suivante : 

Etape n°7 : Identification des vecteurs b, pour (1 <i<M) a I'instants f K+ i en 
appliquant par exemple une technique de separation et d'identification de 
5 sources comme decrit dans les references [2] [3]. 

Etape n°8 : Pour 1<m<M determination des indices i(m) en appliquant 
('equation (4) et en utilisant le vecteur et les vecteurs b, pour (1 <i<M). 

Etape n°9 : Pour 1 <m<M est effectue I'operation a K+ i m= b l(m) 

Etape n»10 : Pour 1<m<Mcalcul du critere U +1 (x,y,z) de facon iterative en 
10 utilisant les equations (14) et (15) et minimisation de U +1 (x,y,z) pour obtenir 
la position (x m ,y m ,z m ) du ni dme emetteur. 

Etape n°11 : Si I'on decide de continuer pour etre plus precis et moins 
ambigue le procede retourne a I'etape n°7. 

Exemple de mise en ceuvre du procede 

15 Les simulations ont ete realisees avec un reseau de N=2 capteurs 

alignes dans I'axe du porteur avec dA=3. Comme d/X=3 une methode 
effectuant une goniometrie aux instants feserait completement ambigue et ne 
permettrait pas de faire par la suite des triangulations pour effectuer la 
localisation de I'emetteur. Sur les figures 5, 6 et 7 correspondant aux criteres 

20 de localisation pour K=3, 7 et 16 est trace le pseudo-spectre L^y) a 
maximiser permettant de determiner la position de I'emetteur dans I'espace 
(x,y). Sachant que si I'emetteur se situe en (xo >yo ) alors /. K (xo,yo)=1 I on en 
deduit que les courbes iso-niveaux Z.K(x,y)=0.99 caracterisent la largeur du 
lobe principale. Remarquant que la precision de localisation depend de la 
25 largeur de ce lobe, on en deduit d'apres les figures 5, 6 et 7 que plus Kest 
important et meilleur sera la precision de localisation. 
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REVINDICATIONS 



1 - Procede de localisation d'une ou de plusieurs sources, la ou lesdites 
sources etant en mouvement par rapport a un reseau de capteurs, le 
procede comportant une etape de separation des sources afin d'identifier les 
vecteurs directeurs associes a la reponse des capteurs a une source 
d'incidence donnee, caracterise en ce qu'if comporte au moins les etapes 
suivantes : 

• associer les vecteurs directeurs a 1m ...a Km obtenus pour le m^me 
emetteur et respectivement pour les instants t r ..t K , 



o localiser le m'^e emetteur a partir des vecteurs 



a im- a Km associes. 



2 - Procede selon la revendication 1 caracterise en ce que I'etape 
□"association comporte au moins les etapes suivantes : 

15 Etape ASE - 1 : initialiser le processus a k=2, 

Etape ASE - 2 : pour 1 <m<M determiner les indices i(m) en utilisant la 
relation d^, b i(m) ) = mn[d(a to a,)], le vecteur a Km et les vecteurs b, 

identifies a I'instant t M pour (1</<M), etablir une fonction p m (t k ) = dCa^, a om ) 

Etape ASE - 3 : pour 1<m<M effectuer Toperation a k+1rri = b i(mh 

20 Etape ASE - 4 : incrementer k<r-k+l et si k<K retourner a I'etape ASE-1 , 

Etape ASE - 5 : a partir de la famille d'instants *={ t,<...< t K }, ainsi 
obtenue, extraire les instants f,- qui n'appartiennent pas a une zone definie 
par la courbe J3 m (t k ) et une zone de tolerance. 
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3 - Procede selon la revendication 1 caracterise en ce que I'etape de 
localisation comporte au moins les §tapes suivantes : 

maximiser un critere de correlation vectorielle normalise L k (x, y, z) dans 
I'espace (x, y, z) de position d'un emetteur avec 

I H I 2 

5 Mx,y,z)= i \ — H — 

Avec 



'\m 



9 Km 



= V K (Xm,ym,Zm) + W K , V^X.y.z) = 



et wk= 



w 



>V 



Km 



oil W/r est le vecteur bruit pourtoutes les positions (x, y, z) d'un emetteur. 
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4 - Procede selon la revendication 3 caracterise en ce que le vecteur b K 
comporte un vecteur representatif du bruit dont les composantes sont 
fonctions des composantes des vecteurs ai m ... aKm- 



15 5 - Procede selon la revendication 3 caracterise en ce qu'il comporte une 
etape ou I'on determine la matrice de covariance R=E[w K w K H ] du vecteur 
bruit et en ce que I'on maximise le critere 

|b/R- l v g (x,j/,z)| 2 

U ' (X,y ' Z) "(b/R- , bJv /c (x^ 5 z) H R-v,(x^ >Z )j 
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6 - Precede selon la revendication 5 caracterise en ce que revaluation du 
critere U{x,y,z) et/ou du critere U'(x,y,z) est recursive. 



7 - Precede selon I'une des revendications 1 a 6 caracterise en ce qu'il 
comporte une etape de comparaison des maximums avec une valeur seuil. 



8 - Precede selon I'une des revendications 1 a 7 caracterise en ce que la 
valeur de Kest initialement fixee a Kq. 



9 - Procede selon I'une des revendications 1 a 8 caracterise en ce que les 
emetteurs a localiser sont mobiles et en ce que le vecteur considere est 
parametre par la position de I'emetteur a localiser et le vecteur vitesse. 
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